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Метод адаптивного контрастирования 
мультиспектральных изображений для 
повышения достоверности нечеткой 
кластеризации 


В статье предложен метод адаптивного повышения контраста мультиспектральных изображений, 
применение которого позволяет повысить достоверность нечеткой кластеризации в случае разбиения 
на малое число классов, характерном при решении задач автоматизированного распознавания символов. 
Представлены экспериментальные результаты применения предложенного метода для обработки КОВ 
снимков железнодорожных вагонов. 


Введение 


В настоящее время одной из актуальных задач является автоматизированное 
распознавание символов, например, номеров железнодорожных вагонов, контейнеров, 
автомобилей и тп. Необходимость таких систем обусловлена сложностью «ручного» 
визуального контроля оператором за процессом перевозки грузов, проездом 
автотранспорта и т.п. в условиях высокой плотности транспортного потока, в то время 
как автоматизированные системы могут легко справиться с этой задачей в режиме 
реального времени. 

Однако существует ряд факторов, усложняющих создание таких систем. Среди них 
можно выделить: возможность нефиксированного положения распознаваемых символов, 
их различный размер и цвет в пределах одного изображения; наличие помех, мешающих 
выделению символов на изображении (например, цемент, просыпанный на номер вагона 
ит.п.), расположение их на плохо освещенной при съемке части объекта. 


Постановка задачи 


Важным этапом в решении задачи распознавания символов является бинаризация 
исходного мультиспектрального изображения. Бинаризация предполагает проведение 
сегментации изображения, в процессе которой его необходимо разделить на 2 кластера, 
к которым будут относиться символы и фон соответственно. Однако в случае 
распознавания номеров железнодорожных вагонов следует учитывать, что кроме 
основного номера на вагоне может располагаться еще и контрольный номер, 
нанесенный существенно меньшим шрифтом и расположенный в его нижней части, 
причем его цвет может отличаться от цвета основного номера. Таким образом, в пределах 
одного изображения символы, требующие распознавания, могут быть изображены 
существенно отличающимися друг от друга цветами. Ситуация может осложняться 
неоднородностью освещения при съемке вагона, что приводит к изменению цвета как в 
пределах одной последовательности символов (номера), так и в пределах одного 
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символа. Следовательно, для повышения чувствительности возникает необходимость 
осуществления кластеризации исходного изображения с разбиением на 3 кластера. 

Применение нечеткой кластеризации непосредственно к исходному изображению 
не позволяет получить удовлетворительных результатов в случае наличия помех и 
низкоконтрастных участков из-за малого числа используемых кластеров, поэтому 
возникает необходимость применения методов повышения контраста и фильтрации. 

Большинство существующих методов повышения контраста и фильтрации [1], [2] 
либо не справляются с поставленной задачей, либо требуют «ручную» подстройку управ- 
ляющих параметров метода под конкретное изображение, что в случае автомати- 
зированной обработки в режиме реального времени не представляется возможным. 

Целью данной статьи является разработка метода адаптивного контрастирования 
мультиспектрального изображения для повышения достоверности его последующей 
нечеткой кластеризации. 


Решение задачи 


Сложность решения поставленной задачи заключается в том, что обработка 
изображений производится автоматизированно, что требует минимизации количества 
управляющих параметров. Кроме того, в пределах одного изображения могут содер- 
жаться фрагменты, очень различные по своим характеристикам: высоко- и низкокон- 
трастные области, которые могут содержать или не содержать помехи. При этом все 
перечисленные особенности изображения используемый метод обработки должен 
учитывать. 

Предложенный метод повышения контраста предполагает адаптивную обработку 
изображения, которое разбивается на неперекрывающиеся окна, и включает в себя 3 
основных этапа, направленные, с одной стороны, на повышение контраста и устранение 
помех на определенных участках исходного изображения, а с другой стороны — на 
сохранение уровня интенсивности высококонтрастных его фрагментов. 

1. Исходное трехдиапазонное изображение, значения цветов которого находятся на 
отрезке [0..1], преобразуется в полутоновое, причем для текущего окна выполняется 
следующее преобразование: 


1 (0.5-йп __ р+тп _р) 


3 
2 (1) 
11 


где СМ у). — пиксель с координатами х, у текущего окна одного из трех диапазонов 


ИСХОДНОГО изображения; У, т пиксель с координатами х,у текущего окна выходного 


полутонового изображения; шп_у -— вектор значений, получаемый на основе 


вычислений среднего значения интенсивности по каждому из трех диапазонов 
текущего окна исходного изображения; ип _р вычисляется по формуле 

. } 1-— ] 

т _р=ши пп (маи) (2) 


У 


где шп_пи, -— среднее значение, вычисляемое для вектора средних значений 


интенсивностей по каждому из трех диапазонов исходного изображения в целом. 
Аналогично вычисляется и значение тп _р: 


ши _р=ши _ УГ), (3) 


где шп _ у, — среднее значение вектора шп _У. 
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2. Производится повышение контраста, полученного на предыдущем этапе полу- 
тонового изображения, состоящее также из трех шагов: 

а) На первом шаге каждое окно входного для этого шага изображения 
подвергается следующему преобразованию: 


м =), (4) 


ху ху 


где \,., — пиксель с координатами х,у текущего окна выходного изображения 


первого шага второго этапа; К — коэффициент, который вычисляется на основе 
статистических характеристик текущего окна входного для этого этапа изображения 
по следующей формуле: 


К =. ны (5) 


где шп _\ - среднее значение интенсивности пикселей текущего окна, входного 
изображения для очередного шага второго этапа, а С вычисляется следующим образом: 


С=(- ша _ №) "9. (6) 


6) В результате выполнения преобразования (4) для каждого пикселя текущего 
окна происходит сдвиг уровня яркости к значениям 1 либо 0 в зависимости от 
коэффициента К, который вычисляется на основе интенсивности пикселя и среднего 
значения интенсивности окна. Однако возможны ситуации, когда это преобразование 
приводит к чрезмерному сдвигу яркости пикселя. Для компенсации возможного 
негативного влияния преобразования (4) и предназначен второй шаг, который 
заключается в применении следующего преобразования к текущему окну изображения: 

Р 
м, =(2,} , (7) 


х,у 


где \ ‚ - пиксель с координатами х,у текущего окна выходного изображения 


второго шага второго этапа; Р - значение, которое вычисляется на основе статисти- 
ческих характеристик текущего окна входного для этого этапа изображения: 


Р=Р", (8) 
где значение Р, вычисляется следующим образом: 
Р =! +$* (пп _ у" —* - 0.5* тах_], (9) 


где С вычисляется по формуле (6), тах_\ — максимальное значение интенсивности 


пикселей текущего окна, входного изображения для этого шага второго этапа, а 5 
вычисляется по формуле 


$ = меп(ти _\№-0.5* тах _\,). (10) 


в) Последний шаг второго этапа позволяет осуществить управляемую коррекцию 
яркости полученного изображения с использованием статистических характеристик 
текущего окна и всего изображения в целом, для осуществления которой текущее 
окно подвергается такому преобразованию: 


шп _\ ) 


Е (о 
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где \*, — пиксель с координатами х,у текущего окна выходного изображения 


последнего шага второго этапа; Р‚ вычисляется по следующей формуле: 
Р, = (1 ша _ на)" +Р_АОЮ, (12) 


где шп_пп - среднее значение яркости по всему полутоновому изображению; 
Р_АОШО - выбираемое эмпирически значение, позволяющее управлять уровнем 


яркости получаемого изображения. 
3. Осуществляется повышение контраста, полученного на предыдущем этапе 
изображения путем преобразований гистограммы для текущего окна: 


(6156; — пой, )* (пах, — пит, ) 


В150 = 


1 


- + т, (13) 
тах, — пт, 

где №15 — вектор значений гистограммы текущего окна изображения, состоящий из 
256 элементов; №15 — новый вектор значений гистограммы обрабатываемого окна; 
тах, пт, — максимальное и минимальное значения интенсивностей, присутствующих 
в обрабатываемом окне соответственно; тах, пи, — новые значения для максимальной 
и минимальной интенсивности окна, которые вычисляются с использованием 


статистических характеристик его гистограммы: 


п, = то ""* (1-4), (14) 


где 4 - доля ненулевых интенсивностей в текущем окне. При этом если выполняется 
условие 0.5 < шш, < 1, то значение шш, корректируется: 


: в я (и, у ) 
шш, =шш, | (15) 
Значение тах, определяется следующим образом: 


| тах —шп _ В 
тах(ти _В,1-ши _В) 


ах =, (16) 
1-4 
где тп _В среднее значение для вектора №15. Если тах, < 0.5, то значение тах, 
корректируется по формуле 
тах, = плах (хи) р (17) 


в противном случае, если выполняется условие шп _В > 0.5, то значение тах, 
изменяется следующим образом: 
тах, = тах ("ами _№), (18) 


Следует отметить, что полученные значения пш, и тах, всегда удовлетворяют 
следующим условиям: 0 < шш, < шш, и тах, < шах, <1. Таким образом, преобразо- 


вание (13) осуществляет масштабированное расширение исходных значений 


гистограммы текущего окна. К полученному вектору №15 применяется следующее 
преобразование: 


в5е = (ве уме (19) 
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. 12 > 
где [15° — вектор значений гистограммы текущего окна выходного изображения; 
БК — коэффициент, вычисляемый на основе статистических характеристик вектора 
№15 по формуле: 
Е 
К — 6156; —шп_ В 
тах(ти _ В — мае _ В) 


(20) 


После этого все интенсивности текущего окна меняются на соответствующие 
им интенсивности, содержащиеся в векторе №15”, в результате чего и формируется 
выходное изображение. 

Экспериментальные результаты были получены на примере обработки 
трехдиапазонных (ВОВ) снимков железнодорожных вагонов размерностью 1600 х 1200 
пикселей, пример которого представлен на рис. 1. На этом изображении основной 
номер изображается черным цветом, а вспомогательный номер, находящийся на низко 
контрастном фрагменте снимка (левый нижний угол) — белым. На рис. 2 представлены 
гистограммы каждой из цветовых составляющих исходного изображения. 


Рисунок 2 — Гистограммы цветовых составляющих исходного КОВ изображения 
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При применении предложенного метода для повышения контраста исходного 
снимка размеры окон для каждого этапа были выбраны 15 х 15, ЗхЗ и 9х9 пикселей, 
соответственно, а параметр Р _ АБП = -0.05. Изображения, полученные после каждого 
этапа обработки исходного снимка, и их гистограммы представлены на рис. 3 — 5 
соответственно. 


Г: г 


ИЗ 
7 ХАЙ 73а 
| ‹ 
=_= | ежззая х10 
В | | 5705358 о 3 


Брута 89 т 


Е АС : 04 05 
а) 6) 
Рисунок 3 — а) полутоновое изображение, полученное после применения 

преобразований 1-го этапа предложенного метода; 6) его гистограмма 
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$%1>< 
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Рисунок 4 — а) изображение, полученное после применения преобразований 
2-го этапа предложенного метода; 6) его гистограмма 
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Нечеткая кластеризация полученного в результате повышения контраста 
изображения производилась методом гибридной кластеризации ЗОМ-ЕСМ, 
описанным в работе [3], результаты которой представлены на рис. 6. Изображение 
разделялось на 3 нечетких кластера, максимальное количество итераций = 2, 
параметры тш=2, пороговое значение ==10“. В процессе кластеризации про- 
изводилось двукратное увеличение числа кластеров перед применением 5ОМ, 
после использования карты Кохонена выбирались 3 кластера с максимальным 
количеством относящихся к ним (имеющим минимальное евклидово расстояние) 
пикселей [4]. 


05 
6) 
Рисунок 5 — а) изображение, полученное после применения преобразований 

3-го этапа предложенного метода; 6) его гистограмма 


Рисунок 6 — Результаты нечеткой кластеризации 
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Применим предложенный метод и в случае обработки КОВ изображений, содер- 
жащих помехи, искажающие распознаваемые символы, то есть так называемые «загря- 
зненные» номера. Пример такого снимка железнодорожного вагона представлен на рис. 7. 
Три символа основного номера на этом изображении частично не видны, остальные же 
символы видны не отчетливо из-за того, что сливаются с фоном, цвет которого искажен 
просыпанным цементом. На рис. 8 представлены результаты нечеткой кластеризации 
изображения, полученного после применения предложенного метода с использованием 
вышеперечисленных управляющих параметров к исходному снимку. Очевидно, что 
уменьшение количества используемых цветов не привело к потере информативности при 
представлении символов номера по сравнению с исходным изображением. 


Рисунок 8 — Результаты нечеткой кластеризации 
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Выводы 


Предложенный адаптивный метод повышения контраста мультиспектральных 
изображений предназначен для получения нового пространства признаков (полуто- 
нового изображения) благодаря применению описанных выше оконных преобразо- 
ваний, направленных на улучшение разделения между яркостями символов и фона, 
что позволяет существенно повысить достоверность нечеткой кластеризации при 
использовании малого числа кластеров. 

К достоинствам предложенного подхода кроме его адаптивности, благодаря 
которой происходит выделение вспомогательных номеров на рассмотренных снимках 
железнодорожных вагонов, можно отнести относительно малое число независящих 
друг от друга управляющих параметров, выбор которых осуществляется при перво- 
начальной настройке на обрабатываемые изображения и не требует дальнейшего 
изменения. 
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Метод адаптивного контрастування мультиспектральних зображень для шдвищення 
достоврност! нечткот кластеризацй 

У статт! запропонований метод адаптивного шдвищення контрасту мультиспектральных зображень, 
застосування якого дозволяе шдвищити вропднсть нечтко! кластеризацй у випадку розбиття на 
мале число клас1в, характерному при розв’язанн! задач автоматизованого розшзнавання символ. 
Представлен! експериментальн! результати застосування запропонованого методу для обробки КОВ 
знмкв зал1зничних вагонйВ. 
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